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[l]  The  probability  distribution  function  (PDF)  of  the  upper 
(0-50  m)  tropical  Pacific  current  speeds  (w),  constructed 
from  hourly  ADCP  data  (1990-2007)  at  six  stations  for 
the  Tropical  Atmosphere  Ocean  project,  satisfies  the  two- 
parameter  Weibull  distribution  reasonably  well  with  different 
characteristics  between  El  Nino  and  La  Nina  events:  In  the 
western  Pacific,  the  PDF  of  w  has  a  larger  peakedness  during 
the  La  Nina  events  than  during  the  El  Nino  events;  and  vice 
versa  in  the  eastern  Pacific.  Flowever,  the  PDF  of  w  for  the 
lower  layer  (100-200  m)  does  not  fit  the  Weibull  distribution 
so  well  as  the  upper  layer.  This  is  due  to  the  different 
stochastic  differential  equations  between  upper  and  lower 
layers  in  the  tropical  Pacific.  For  the  upper  layer,  the 
stochastic  differential  equations,  established  on  the  base  of 
the  Ekman  dynamics,  have  analytical  solution,  i.e.,  the 
Rayleigh  distribution  (simplest  form  of  the  Weibull 
distribution),  for  constant  eddy  viscosity  K.  Knowledge  on 
PDF  of  w  during  the  El  Nino  and  La  Nina  events  will  improve 
the  ensemble  horizontal  flux  calculation,  which  contributes 
to  the  climate  studies.  Citation:  Chu,  P.  C.  (2008),  Probability 
distribution  function  of  the  upper  equatorial  Pacific  current  speeds, 
Geophys.  Res.  Lett,  35,  L12606,  doi:10.1029/2008GL033669. 


[4]  The  w-PDF  has  been  investigated  thoroughly  in  the 
atmosphere.  For  example,  the  vc-PDF  for  the  surface  winds 
is  well  represented  by  the  two-parameter  Weibull  distribu¬ 
tion  [e.g.,  Monahan,  2006],  However,  the  w-PDF  has  not 
been  investigated  in  the  oceans.  To  fill  this  gap,  we  use 
hourly  Acoustic  Doppler  Current  Profiler  (ADCP)  data 
(1990-2007)  at  all  the  six  stations  during  the  Tropical 
Atmosphere  Ocean  (TAO)  project  [MePhaden  et  al., 
1998]  in  this  study  to  construct  the  observational  w-PDF 
for  various  layers.  The  purpose  is  to  identify  which  theo¬ 
retical  PDF  should  be  used  for  the  current  speed  for  general 
and/or  special  events  such  as  El  Nino  and  La  Nina.  Special 
characteristics  of  the  statistical  parameters  such  as  mean, 
standard  deviation,  skewness,  and  kurtosis  will  also  be 
identified.  The  rest  of  the  paper  is  organized  as  follows. 
Section  2  describes  the  data.  Section  3  presents  theoretical 
background  for  determining  the  w-PDF  for  upper  oceans. 
Section  4  shows  the  basic  characteristics  of  the  four 
parameters  of  the  Weibull  distribution.  Section  5  constructs 
the  observational  w-PDF  from  the  TAO-ADCP  data  and 
discusses  the  basic  features  of  the  four  parameters.  Section  6 
presents  the  conclusions. 


1.  Introduction 


2.  Data 


[2]  Tropical  Pacific  Ocean  contributes  significantly  to  the 
global  redistribution  of  heat  necessary  to  maintain  the  earth’s 
thermal  equilibrium.  For  example,  the  El  Nino  (La  Nina) 
phenomenon,  an  eastward  (westward)  shift  of  warm  (cool) 
water,  is  a  key  component  of  global  interannual  climate 
variability.  In  connection  with  the  El  Nino  and  La  Nina 
phenomena,  impact  of  upper  tropical  Pacific  on  global  climate 
changes  has  attracted  large  attention.  During  the  El  Nino  or 
La  Nina  event,  equatorial  currents  take  an  active  role  in 
redistribution  of  heat  that  changes  the  sea  surface  temperature 
(SST)  and  in  turn  affects  the  atmospheric  general  circulations. 

[3]  Surface  layer  horizontal  fluxes  of  momentum,  heat, 
water  mass  and  chemical  constituents  are  typically  nonlin¬ 
ear  in  the  speed  [ Lozano  et  al,  1996;  Galanis  et  al.,  2005], 
so  the  space  or  time  average  flux  is  not  generally  equal  to 
the  flux  that  would  be  diagnosed  from  the  averaged  current 
speed.  In  fact,  the  average  flux  will  generally  depend  on 
higher-order  moments  of  the  current  speed,  such  as  the 
standard  deviation,  skewness,  and  kurtosis.  From  both 
diagnostic  and  modeling  perspectives,  there  is  a  need  for 
parameterizations  of  the  probability  distribution  function 
(PDF)  of  the  current  speed  w  (called  w-PDF  here). 
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[5]  The  upper  layer  ocean  current  dataset  consists  of 
hourly  ADCP  moorings  from  the  TAO  project  at  the  six 
stations  (147°E,  156°E,  165°E,  170°W,  140°W,  110°W) 
along  the  equator  with  the  vertical  resolution  of  5  m. 
Temporal  data  coverage  varies  from  station  to  station: 
1993-2004  with  interruptions  during  1997-1998  and 
2001-2003  on  147°E,  1990-1995  on  156°E,  1990-2007 
with  interruption  during  1993-1994  on  165°E,  1999-2006 
with  interruption  1992-1993  on  170°W,  1995-2007  on 
140°W,  and  1995-2007  with  interruption  in  2004  on 
110°W.  The  hourly  horizontal  velocity  data  (u,  v )  can  be 
downloaded  from  the  website:  http://www.pmel.noaa.gov/ 
tao /.  The  current  speed  (w)  is  calculated  from  ( u ,  v). 


3.  Theoretical  Background 


[6]  Let  (x,  y)  be  the  horizontal  coordinates  and  z  be  the 
vertical  coordinate.  The  corresponding  horizontal  velocity 
components  are  represented  by  (it,  v).  Large-scale  horizontal 
momentum  equation  is  written  by 


du 

dt 


1  dx 

p  dz  ’ 


(1) 


dv 

dt 


+f(u-ug) 


1  dY 
p  dz  ’ 


(2) 
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where  ( X ,  Y)  are  the  horizontal  stresses;  p  is  the  pressure;  p 
is  the  density;  /  is  the  Coriolis  parameter;  and  ( ug ,  vg)  are 
geostrophic  velocity  components  defined  by 


1  dp  1  dp 

lg  =  ~fp  dy ’  Vg=Tp~d~x' 


(3) 


[9]  For  the  sake  of  convenience,  we  assume  that  the 
residual  between  the  Ekman  transport  (UE,  Ve)  and  surface 
wind  stress  does  not  depend  on  the  horizontal  current  vector 
(w,  v).  Away  from  the  equator,  this  approximation  is  similar 
to  a  small  Rossby  number  approximation  [Gz7/,  1982].  If  the 
forcing  (A„,  Av)  is  fluctuating  around  some  mean  value, 


[7]  Integrating  ( 1 )  and  (2)  from  the  ocean  surface  (z  =  0) 
to  a  constant  scale  depth  (h)  of  surface  mixed  layer  leads  to 

d^-JVE  =  -p(X0-X_h),  (4) 

%+JUb=-{Yo-Y-M),  (5) 

at  p 

where 

(U,V)=  f  (u,v)dz  =  (hu,hv),  (6) 


,° 

{UE,VE)  =  (u  -  ug,v  -  vg)dz,  (7) 

J-h 

with  ( u ,  v)  the  vertical  means  of  ( u ,  v),  ( X0 ,  Y0)  =  (tx,  ry)  the 
surface  wind  stress  components;  and  (X_h,  Y_h)  the  stress 
components  at  z  =  —  h,  which  is  calculated  by 


A„W  =  <A„}+  Wi(t)h£,  \v{t)  =  <AV>  +  Wi{t)KZ,  (12) 

where  the  angle  brackets  represent  ensemble  mean  and  the 
fluctuations  are  taken  to  be  isotropic  and  white  in  time: 

(Wt{h)Wj{t2))  mSySih  -t2),  (13) 

with  a  strength  that  is  represented  by  E.  Note  that  the 
Ekman  transport  is  detennined  by  the  surface  wind  stress 
for  time-independent  case,  and  therefore  the  ensemble  mean 
values  of  (A„,  Av)  are  zero, 

<A„)=0,  (Av)  =  0.  (14) 

[10]  Substitution  of  (12)— (14)  into  (9)  and  (10)  gives 

du  K  ■ 

+  (15) 

|  — p'+«')£,  (16) 


A_,,  K  Y_h  K 

- ~~TU’  - ~-Tv. 

p  h  p  h 


(8) 


[8]  Here,  we  assume  that  the  horizontal  velocity  is  much 
weaker  below  the  mixed  layer  than  in  the  mixed  layer.  K  is 
the  eddy  viscosity.  Substitution  of  (6)-(8)  into  (4)  and  (5) 
gives 


du  l  K 
~di~hAu~¥u’ 


(9) 


dv  1  K 

~di~hAv~JYu’ 


(10) 


which  is  a  set  of  stochastic  differential  equations  for  the 
surface  current  vector.  The  joint  PDF  of  (m,  v)  satisfies  the 
Fokker-Planck  equation, 


dp  _  (d1p 

dt~  \  2)  \du^ 


(17) 


which  is  a  linear  second-order  partial  differential  equation 
with  the  depth  scale  ( h )  taken  as  a  constant.  Transforming 
from  the  orthogonal  coordinates  (m,  v)  to  the  polar 
coordinates  (w,  (p)  respectively  the  current  speed  and 
direction, 


where 


u  =  w  cos  <p,  v  =  w  sin  (p. 


(18) 


A u  =jVE  H ,  Av  =  — fUE  +  — 

P  P 


(11)  [n]  The  joint  PDF  of  (u,  v)  is  transformed  into  the  joint 

PDF  of  (w,  (p ), 


represent  the  residual  between  the  Ekman  transport  and 
surface  wind  stress.  With  absence  of  horizontal  pressure 
gradient,  e.g.,  ug  =  vg  =  0,  Equations  (9)  and  (10)  reduce  to 
the  commonly  used  wind-forced  slab  model  [e.g.,  Pollard 
and  Millard,  1970], 


p(u,v)dudv  =  p(u,v)wdwd(p  =  p(w,(p)dwd(p.  (19) 

[12]  Integration  of  (19)  over  the  angle  <p  from  0  to  27r 
yields  the  marginal  PDF  for  the  current  speed  alone, 


du 

di 


K 

¥-u 


p2tt 

p{w)=  /  p(w,  tp)dtp.  (20) 

Jo 


dv 

dt 


—fit  T 


ph 


K 


[13]  For  a  constant  eddy  viscosity  ( K)  at  z  =  —h,  the 
steady  state  solution  of  equation  (17)  is  given  by 
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p{u,  v)  =  A  exp 


K 

Wh2 


(i U 2  +  V2) 


(21) 


where  A  is  a  normalization  constant.  Substitution  of  (21) 
into  (19)  and  use  of  (20)  yield 


p(w)  =  2ivAwexpl  — 


Kw2 

Wh2 


(22) 


with 


where  T  is  the  gamma  function.  The  parameters  a  and  b  can 
be  inverted  [Monahan,  2006]  from  (26)  and  (27), 


mean(w) 

std(w) 


mean(w) 

r(i +  i/&)‘ 


(30) 


[i6]  The  skewness  and  kurtosis  depend  on  the  parameter 
b  only  [see  (26)  and  (27)]  for  the  Weibull  distribution.  The 
relationship  between  the  kurtosis  and  skewness  can  be 
detennined  from  (28)  and  (29). 


J  p(w)dw  =  1.  (23) 

o 

[14]  Substitution  of  (22)  into  (23)  leads  to  the  Rayleigh 
distribution 


p{w) 


2  w 
2  exp 
a 2 


_  T,h 


(24) 


with  the  scale  parameter  a.  The  basic  postulation  of  constant 
K  may  not  be  met  always  at  the  upper  ocean.  Hence  we 
require  a  model  that  can  meet  the  twin  objectives  of  (a) 
accommodating  Rayleigh  distribution  whenever  the  basic 
hypothesis  (constant  K)  that  justifies  it  is  satisfied  and  (b) 
fitting  data  under  more  general  conditions.  This  requirement 
is  supposed  to  be  satisfied  by  the  Weibull  probability 
density  function, 


p(w)  = 


/WN*-1 

(w\~ 

-)  exp 

-(-) 

\aJ 

\aJ 

(25) 


where  the  parameters  a  and  b  denote  the  scale  and  shape  of 
the  distribution.  This  distribution  has  been  recently  used  in 
investigating  the  ocean  model  predictability  by  Ivanov  and 
Chu  [2007a,  2007b], 


4.  Parameters  of  Weibull  Distribution 

[15]  The  four  parameters  (mean,  standard  deviation, 
skewness,  and  kurtosis)  of  the  Weibidl  distribution  are 
calculated  by  [Johnson  et  al.,  1994] 

mean(w)  aT  M  +  — j  ,  (26) 


std(w) 


-T2 


(27) 


5.  Observational  w-PDFs 

[17]  The  data  depicted  in  Section  2  are  used  to  construct 
observational  w-PDF  (i.e.,  histograms)  for  the  two  stations 
(165°E,  1 10°W)  along  the  equator  in  upper  oceans  (0-50  m) 
for  the  whole  period  (1990-2007),  major  El  Nino  events 
(May  91-Jul  92,  Dec  92-Jun  93,  Jul  94-Mar  95,  May  97- 
Apr  98,  May  02-Mar  03,  Jul  04-Feb  05,  Sep  06-Jan  07), 
and  major  La  Nina  events  (Sep  95 -Mar  96,  Jul  98-Jun  00, 
Oct  00-Feb  01,  Aug-Dec  07).  It  is  found  that  the  observa¬ 
tional  w-PDF  fit  the  two-parameter  Weibull  distribution 
reasonably  well  in  all  occasions  (Figure  1).  At  165°E 
(western  Pacific)  the  w-PDF  has  a  largest  peakedness  with 
a  lower  mode  (0.25  m  s-1)  during  the  La  Nina  events,  a 
medium  peakedness  for  the  whole  period,  and  smallest 
peakedness  with  a  higher  mode  (0.4  m  s-1)  during  the 
El  Nino  events;  and  vice  versa  at  110°W  (eastern  Pacific), 
the  w-PDF  has  a  largest  peakedness  with  a  lower  mode 
(0.4  m  s~  1 )  during  the  El  Nino  events,  a  medium  peaked¬ 
ness  for  the  whole  period,  and  smallest  peakedness  wit  a 
higher  mode  (0.65  m  s-1)  during  the  La  Nina  events. 
From  the  western  to  eastern  Pacific,  the  mode  is  compa¬ 
rable  (0.4  m  s-1)  during  the  El  Nino  events,  and  increases 
from  0.25  m  s'1  to  0.65  m  s-1  during  the  La  Nina  event. 
Such  flip  pattern  of  w-PDF  between  the  eastern  and 
western  Pacific  during  El  Nino/La  Nina  events  may  imply 
the  importance  of  equatorial  current  systems  in  the  El 
Nino-Southern  Oscillation  phenomenon. 

[is]  The  four  parameters  (mean,  standard  deviation, 
skewness,  and  kurtosis)  can  also  be  calculated  from  the 
observational  data  (w), 


1  N  / - - 

mean(w)  =  —  ^  w,-,  std(w)  =  y  mean[w  —  mean(w)]~, 


skew(w)  = 


mean<  \w  —  mean w) 


’} 


std  (w) 


kurt(w)  =  ■ 


mean!  w  —  mean!  wl 


std4(w) 


--3, 


(31) 


skew(w) 


r(i+|)-3r(i+i)r(i+|)  +  2r3(1+i) 
[r(i  +  i)-r2(i+i)]3/2 

(28) 


for  any  location  (station,  depth)  where  the  ADCP  measure¬ 
ments  were  taken.  Thus,  a  four-parameter  dataset  has  been 
established  each  location.  The  scatter  diagrams  are  drawn 
for  all  six  stations  (147°E,  156°E,  165°E,  170°W,  140°W, 


kurt(w) 


r(i+|)^4r(i+i)r(i+|)  +  6r2(i+i)r(i  +  |)-3r4(i+i) 

[r(i  +  i)-r2(i  +  i)]2 


(29) 
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Figure  1.  Comparison  between  observational  u-PDFs  (i.e.,  histogram,  dashed  curve)  constructed  from  the  TAO-ADCP 
data  and  Weibull  distributions  (solid  curve)  for  upper  layer  (0-50  m)  along  the  equator  at:  (left)  165°E,  and  (right)  110°W 
with  the  upper  panels  for  the  whole  period,  the  middle  panels  for  the  major  El  Nino  events,  and  the  lower  panels  for  the 
major  La  Nina  events. 


and  110°W)  along  the  equator  in  upper  oceans  (0-50  m) 
and  sub-layer  (100-200  m)  for  the  whole  period  (1990- 
2007),  and  El  Nino/La  Nina  events.  We  may  use  the 
relationships  between  the  skewness  and  the  mean/std  ratio 
(representing  the  parameter  b.  Figure  2)  and  between  the 
kurtosis  and  the  skewness  (Figure  3)  to  identify  the  fitness 
of  the  Weibull  distribution  for  observational  w-PDFs  for  the 
whole  period  and  the  El  Nino/La  Nina  events.  The  solid 
curve  on  Figures  1-3  shows  the  relationship  for  a  Weibull 
variable. 

[19]  For  the  observational  n'-PDF  in  the  upper  layer 
(0-50  m),  the  skcw(vv)  is  evidently  a  concave  function  of 
the  ratio  mean(w)/std(w)  (the  same  as  the  Weibull  distri¬ 
bution),  such  that  the  theoretical  function  is  positive  for 
small  values  of  this  ratio  and  negative  for  large  values. 
Flowever,  the  ratio  mean(w)/std(w)  is  always  less  than  2  for 
the  whole  period  and  El  Nino  events,  and  less  than  2.5  for 
the  La  Nina  events.  The  skewness  is  generally  positive  in 


the  upper  tropical  Pacific  for  all  occasions  (Figure  2,  top). 
Similarly,  the  relationship  between  skew(w)  and  kurt(w)  in 
the  observations  is  also  similar  to  that  for  a  Weibull 
variable  (Figure  3).  The  agreement  between  the  moment 
relationships  in  the  upper  layer  (0-50  m)  TAO-ADCP  data 
and  those  for  a  Weibull  variable  reinforces  the  conclusion 
that  these  data  are  Weibull  to  a  good  approximation,  which 
is  not  affected  by  the  large-scale  processes  such  as  El  Nino 
and  La  Nina  events. 

[20]  The  scatter  diagrams  of  skew(w)  versus  ratio 
mean(w)/std(w)  and  kurt(w)  versus  skew(w)  for  the  sub¬ 
layer  (100-200  m)  (bottom  plots  of  Figures  2  and  3)  show 
that  the  w-PDF  for  the  lower  layer  (100-200  m)  does  not  fit 
the  Weibull  distribution  so  well  as  the  upper  layer.  Evident 
difference  is  found  between  the  observational  vv-PDF  and 
the  Weibull  distribution.  Such  difference  can  be  explained 
as  follows.  For  the  upper  layer,  the  horizontal  velocity 
satisfies  (9)  and  (10).  This  leads  to  that  the  current  speed 
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Figure  2.  Relationship  between  the  ratio  mean(w)/std(w)  and  skew(w)  for  the  observational  w-PDFs  (dots)  from 
TAO-ADCP  at  the  six  stations  and  the  Weibull  distribution  (solid  curve)  during  (left)  the  whole  period,  (middle) 
major  El  Nino  events,  and  (right)  major  La  Nina  events  with  the  upper  panels  for  the  upper  layer  (0-50  m),  and  the 
lower  panels  for  the  sub-layer  (100-200  m). 
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Figure  3.  Relationship  between  kurt(w)  and  skcw(vv)  for  the  observational  vv-PDFs  (dots)  from  TAO-ADCP  at  the  six 
stations  and  the  Weibull  distribution  (solid  curve)  during  (left)  the  whole  period,  (middle)  major  El  Nino  events,  and 
(right)  major  La  Nina  events  with  the  upper  panels  for  the  upper  layer  (0-50  m),  and  the  lower  panels  for  the  sub-layer 
(100-200  m). 
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(w)  satisfies  the  Rayleigh  distribution  for  the  constant  eddy 
viscosity  K  and  may  extend  to  more  general  Weibull 
distribution  for  non-constant  K.  For  the  lower  layer,  the 
horizontal  velocity  does  not  satisfy  (9)  and  (10).  This  may 
cause  the  deviation  of  w-PDF  from  the  Weibull  distribution. 

6.  Conclusions 

[21]  This  study  has  investigated  the  probability  distribu¬ 
tion  function  of  the  current  speeds  (w)  in  upper  tropical 
Pacific  observationally,  using  long-term  (1990-2007) 
hourly  ADCP  data  at  six  stations  during  the  TAO  project; 
and  theoretically,  using  a  stochastic  model  derived  using 
boundary  layer  physics.  The  following  results  were 
obtained. 

[22]  (1)  Probability  distribution  function  of  the  current 
speeds  (w)  satisfies  the  two-parameter  Weibull  distribution 
in  the  upper  tropical  Pacific  Ocean  (0-50  m)  and  does  not 
satisfy  the  two-parameter  Weibull  distribution  in  lower 
tropical  Pacific  (100-200  m).  In  the  upper  tropical  Pacific 
with  a  constant  eddy  viscosity  K ,  the  probability  distribution 
function  satisfies  a  linear  second-order  partial  differential 
equation  (i.e.,  the  Fokker-Planck  equation)  with  an  analyt¬ 
ical  solution  -  the  Rayleigh  distribution  (special  case  of  the 
2  parameter  Weibull  distribution).  The  stochastic  differential 
equations  are  different  between  upper  and  lower  layers  in 
tropical  Pacific  (i.e.,  making  the  corresponding  Fokker- 
Planck  equation  different),  which  causes  the  probability 
distribution  function  different. 

[23]  (2)  The  w-PDF  in  the  upper  ocean  (0-50  m)  satisfies 
the  two-parameter  Weibull  distribution  reasonably  well  all 
the  time  with  different  characteristics  between  El  Nino  and 
La  Nina  events.  In  the  western  Pacific,  the  w-PDF  has  a 
larger  peakedness  during  the  La  Nina  events  than  during  the 
El  Nino  events;  and  vice  versa  in  the  eastern  Pacific.  From 
the  western  to  eastern  Pacific,  the  mode  is  comparable 
(0.4  m  s^1)  during  the  El  Nino  events,  and  increases  from 
0.25  m  s  1  to  0.65  m  s'  1  during  the  La  Nina  event.  Such 
flip  pattern  of  w-PDF  between  the  eastern  and  western 
Pacific  during  El  Nino/La  Nina  events  may  imply  the 
importance  of  equatorial  current  systems  in  the  El  Nino- 
Southem  Oscillation  phenomenon. 

[24]  (3)  Four  moments  of  w  (mean,  standard  deviation, 
skewness,  kurtosis)  have  been  characterized.  It  was  found 
that  the  relationships  between  mean(w)/std(w)  and  skew(w) 
and  between  skew(w)  and  kurt(w)  from  the  data  are  in  fairly 
well  agreement  with  the  theoretical  Weibull  distribution  for 
the  upper  (0-50  m)  tropical  Pacific  for  the  whole  period  as 
well  as  El  Nino/La  Nina  events,  but  are  not  in  well 


agreement  with  the  theoretical  Weibull  distribution  for  the 
lower  (below  100  m  depth)  tropical  Pacific.  The  ADCP  data 
also  show  that  the  ratio  mean(w)/std(w)  is  generally  less 
than  2  for  the  whole  period  and  El  Nino  events,  and  less 
than  2.5  for  the  La  Nina  events.  The  skewness  is  generally 
positive  in  the  upper  (0-50  m)  tropical  Pacific. 

[25]  (4)  A  primary  motivation  for  the  study  of  the 
probability  distribution  of  upper  layer  ocean  current  speeds 
is  the  role  these  distributions  play  in  the  computation  of 
spatially  and/or  temporally  averaged  horizontal  fluxes  of 
momentum,  heat,  water  mass  and  chemical  constituents. 
The  Weibull  distribution  provides  a  good  empirical  approx¬ 
imation  to  the  PDF  of  w,  which  is  not  affected  by  the  large- 
scale  processes  such  as  El  Nino  and  La  Nina  events.  This 
presents  the  possibility  of  improving  the  representation  of 
the  horizontal  fluxes  that  are  at  the  heart  of  the  coupled 
physical-biogeochemical  dynamics  of  the  marine  system. 
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